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1  Objectives 


•  Klein  tunneling,  abnomal  electron  paths,  and  extreme  conductance  fluctuations  in  graphene  quantum 
point  contacts. 

•  Control  of  quantum  transport  in  nanostructures  by  classical  chaos. 

•  Relativistic  quantum  Darwinism  and  applications. 

The  AFOSR  support  helped  create  a  new  field  of  interdisciplinary  research:  Relativistic  Quantum 
Chaos,  which  studies  the  relativistic  quantum  manifestations  of  classical  chaos  with  applications  in  graphene 
and  alternative  two-dimensional  Dirac  material  systems. 


2  List  of  Publications 

1.  X.  Ni,  L.  Huang,  Y.-C.  Lai,  and  L.  M.  Pecora,  “Effect  of  chaos  on  relativistic  quantum  tunneling,” 
Europhysics  Letters  98,  50007,  1-6  (2012).  This  paper  was  highlighted  in  Europhysics  News,  Vol.  43, 
No.  6  (2012). 

2.  L.  Ying,  L.  Huang,  Y.-C.  Lai,  and  C.  Grebogi,  “Conductance  fluctuations  in  graphene  systems:  The 
relevance  of  classical  dynamics,”  Physical  Review  B  85,  245448,  1-7  (2012). 

3.  X.  Ni,  L.  Huang,  Y.-C.  Lai,  and  C.  Grebogi,  “Scarring  of  Dirac  fermions  in  chaotic  billiards,”  Physical 
Review  E  86,  015702,  1-8  (2012). 

4.  X.  Ni,  Y.-C.  Lai,  and  W.-X.  Wang,  “Emergence  of  scaling  associated  with  complex  branched  wave 
structures  in  optical  medium,”  Chaos  22,  043116,  1-11  (2012). 

5.  H.-Y.  Xu,  L.  Huang,  Y.-C.  Lai,  and  C.  Grebogi,  “Cbiral  scars  in  chaotic  Dirac  fermion  systems,” 
Physical  Review  Letters  110,  064102,  1-5  (2013).  This  work  was  selected  by  PRL  Editors  as  an 
Editors’  Suggestion. 

6.  L.  Huang,  R.  Yang,  Y.-C.  Lai,  and  D.  K.  Eerry,  “Lead-position  dependent  regular  oscillations  and  ran¬ 
dom  fluctuations  of  conductance  in  graphene  quantum  dots,”  Journal  of  Physics:  Condensed  Matters 
25,  085502,  1-7  (2013). 

7.  L.  Ying,  L.  Huang,  Y.-C.  Lai,  and  Y.  Zhang,  “Effect  of  geometrical  rotation  on  conductance  fluctua¬ 
tions  in  graphene  quantum  dots,”  Journal  of  Physics:  Condensed  Matters  25,  105802,  1-7  (2013). 

8.  R.  Yang,  L.  Huang,  Y.-C.  Lai,  C.  Grebogi,  and  L.  M.  Pecora,  “Harnessing  quantum  transport  by 
transient  chaos,”  Chaos  23,  013125,  1-9  (2013). 

9.  G.-L.  Wang,  L.  Ying,  Y.-C.  Lai,  and  C.  Grebogi,  “Quantum  chaotic  scattering  in  graphene  systems  in 
the  absence  of  invariant  classical  dynamics,”  Physical  Review  E  87,  052908,  1-9  (2013). 

10.  X.  Ni,  L.  Huang,  L.  Ying,  and  Y.-C.  Lai,  “Relativistic  quantum  tunneling  of  Dirac  fermions  in  nonhy- 
perbolic  chaotic  systems,”  Physical  Review  B  87,  224304,  1-10  (2013). 

11.  G.-L.  Wang,  L.  Huang,  Y.-C.  Lai,  and  C.  Grebogi,  “Nonlinear  dynamics  and  quantum  entanglement 
in  optomechanical  systems,”  Physical  Review  Letters  112,  110406,  1-6  (2014). 

12.  L.  Ying,  Y.-C.  Lai,  and  C.  Grebogi,  “Quantum  manifestation  of  synchronization  transition  in  optome¬ 
chanical  systems,”  Physical  Review  A  90,  053810,  1-5  (2014). 

13.  L.  Ying,  G.-L.  Wang,  L.  Huang,  and  Y.-C.  Lai,  “Quantum  chaotic  tunneling  in  graphene  systems  with 
electron-electron  interactions,”  Physical  Review  B  90,  224301,  1-14  (2014). 
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14.  L.  Huang,  Y.-C.  Lai,  H.-G.  Luo,  and  C.  Grebogi,  “Universal  formalism  of  Fano  resonance,”  AIP 
Advances  5,  017137,  1-18  (2015). 

15.  Y.-Y.  Park,  Y.-H.  Do,  S.  Altmeyer,  Y.-C.  Lai,  and  G.-W.  Lee,  “Early  effect  in  time-dependent,  high¬ 
dimensional  nonlinear  dynamical  systems  with  multiple  resonances,”  Physical  Review  E  91,  022906, 
1-7  (2015). 

16.  H.-Y.  Xu,  L.  Huang,  Y.-C.  Lai,  and  C.  Grebogi,  “Superpersistent  currents  and  whispering  gallery 
modes  in  relativistic  quantum  chaotic  systems,”  Nature  Scientific  Reports  5,  8963,  1-8  (2015). 

17.  S.  Altmeyer,  Y.-H.  Do,  and  Y.-C.  Lai,  “Transition  to  turbulence  in  Taylor-Couette  ferrofluidic  flow,” 
Nature  Scientific  Reports  5,  10781,  1-12  (2015). 

18.  S.  Altmeyer,  Y.-H.  Do,  and  Y.-C.  Lai,  “Ring  bursting  behavior  en  route  to  turbulence  in  quasi  two- 
dimensional  Taylor-Couette  flows,”  submitted  to  Physical  Review  E. 

19.  P.-P.  Li,  Z.-G.  Huang,  L.  Huang,  and  Y.-C.  Lai,  “Lyapunov  exponent  of  Sinai  billiard,”  submitted  to 
EuroPhysics  Eetters. 

20.  R.  Bao,  L.  Huang,  Y.-C.  Lai,  and  C.  Grebogi,  “Conductance  fluctuations  in  chaotic  bilayer  graphene 
quantum  dots,”  submitted  to  Physical  Review  E. 


3  Accomplishments  and  New  Findings 

3.1  Solutions  of  Dirac  equation,  relativistic  quantum  tunneling  and  scarring 

3.1.1  General  method  for  solving  the  Dirac  equation  in  closed  domains  and  relativistic  quantum 
scarring 

Given  a  closed  Hamiltonian  system  that  exhibits  fully  developed  chaos  in  the  classical  limit,  one  might 
expect  the  quantum  wavefunctions  associated  with  various  eigenstates  to  be  more  or  less  uniform  in  the 
physical  space.  This  had  been  thought  to  be  quite  natural  because  of  ergodicity  associated  with  chaos  in 
the  classical  phase  space.  The  notion  of  uniform  wavefunctions  was  nevertheless  proven  to  be  wrong  about 
three  decades  ago,  where  strongly  non-uniform  eigenfunctions  were  discovered  by  McDonald  and  Kaufman 
in  their  study  of  the  eigen- solutions  of  the  Schrodinger  equation  in  the  chaotic  stadium  billiard.  A  systematic 
study  was  subsequently  carried  out  by  Heller,  who  established  the  striking  tendency  for  wavefunctions  to 
concentrate  about  classical  unstable  periodic  orbits,  which  he  named  quantum  scars.  Semiclassical  theory 
was  then  developed  by  Bogomolny  and  Berry,  providing  a  general  understanding  of  the  physical  mechanism 
of  quantum  scars.  It  should  be  noted  that,  the  phenomenon  was  deemed  counterintuitive  and  surprising 
solely  because  of  chaos,  as  the  phase  space  of  an  integrable  system  is  not  ergodic  so  that  the  quantum 
wavefunctions  are  generally  not  expected  to  be  uniform.  Quantum  scarring  is  one  of  the  most  remarkable 
phenomena  in  modern  physics,  and  has  been  an  active  area  of  research. 

Previous  works  on  scarring  focused  on  non-relativistic  quantum  systems  described  by  the  Schrodinger 
equation.  A  fundamental  issue  was  whether,  in  a  closed  chaotic  domain,  relativistic  quantum  particles 
obeying  the  Dirac  equation  can  scar.  This  issue  was  partially  addressed  by  Pi’s  group  in  the  context  of 
closed  graphene  confinements,  where  signatures  of  quantum  scars  were  identified  in  the  patterns  of  the  local 
density  of  states.  The  framework  used  in  the  study,  however,  was  based  on  the  tight-binding  Hamiltonian  and 
the  non-equilibrium  Green’s  function  derived  still  from  the  Schrodinger  equation.  The  question  of  whether 
truly  Dirac  fermions,  particles  strictly  obeying  the  Dirac  equation,  can  scar  in  chaotic  billiards  had  not  been 
addressed. 
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Figure  1 :  Examples  of  eigenstates  of  the 
Dirac  equation  for  a  bow-tie  chaotic  bil¬ 
liard.  The  maximum  height  (vertical  dis¬ 
tance  from  tip  to  base)  is  set  to  1 ,  and  the 
distance  between  two  tips  is  2.  Panels  (a- 
d)  are  for£  =  36.1335,  43.0190,  27.9300, 
and  36.2729,  respectively.  Panels  (a)  and 
(b)  show  the  first  component  of  the  Dirac 
spinor,  while  (c)  and  (d)  show  the  second 
component. 


(a)  , 


(b) 


(c) 


(d) 


Figure  2:  Examples  of  eigenstates  of  the 
Dirac  equation  for  the  chaotic  “Africa” 
billiard.  The  domain  is  confined  within  a 
rectangular  box  x  S  [—0.99, 1.35]  and  y  G 
[—1.22,0.92].  Panels  (a-d)  are  for  E  = 
15.0468,  32.7638,  13.4838,  and  15.2402, 
respectively.  Panels  (a)  and  (b)  show  the 
first  spinor  component,  while  (c)  and  (d) 
show  the  second  spinor  component. 


Prior  to  our  work,  a  general  method  for  solving  completely  the  Dirac  equation  in  closed  system  of 
arbitrary  geometry  had  not  existed.  The  main  reason  was  that,  although  the  Dirac  equation  is  the  cornerstone 
of  relativistic  quantum  mechanics  and  quantum  electrodynamics,  solutions  were  focused  on  free  space  and 
perturbation  types  in  situations  where  relativistic  quantum  behaviors  occur.  Prior  to  the  discoveries  of 
graphene  and  topological  insulators,  it  was  not  generally  thought  that  relativistic  quantum  mechanics  would 
be  relevant  to  solid  state  devices.  As  a  result,  there  had  been  little  interest  in  studying  the  Dirac  equation 
in  finite  domains.  The  only  exception  was  the  work  of  Berry  and  Mondragon,  who  in  1987  developed  a 
boundary-integral  type  of  method  to  solve  the  energy  levels  (eigenvalues)  of  a  chaotic  neutrino  billiard.  To 
obtain  the  eigenfunctions,  a  closed  form  of  the  boundary  of  the  domain  is  needed. 

Supported  by  AFOSR,  we  developed  a  general  and  efficient  method  to  solve  the  Dirac  equation  for 
massless  fermions  in  two-dimensional  closed  domains.  An  obstacle  to  obtaining  a  complete  solution  of 
the  Dirac  equation,  which  includes  both  eigenvalues  and  eigenfunctions,  was  the  proper  handling  of  the 
boundary  conditions.  We  articulated  an  efficient  discretization  scheme  and  a  physically  meaningful  approach 
to  treating  the  boundary  conditions,  converting  the  Dirac  equation  into  a  set  of  matrix  equations.  In  our 
method,  the  physical  symmetries  of  the  system  are  well  preserved.  To  validate  our  method,  we  considered 
three  types  of  representative  geometric  confinements,  which  include  domains  that  generate  integrable  or 
chaotic  motions  in  the  classical  limit,  and  calculated  the  complete  spectra  of  eigenvalues  and  the  associated 
eigenvector  sets.  In  particular,  in  the  case  of  integrable  geometries  for  which  analytic  predictions  of  the 
eigenvalues  and  eigenvectors  are  available,  we  obtained  excellent  agreement  between  the  numerical  and 
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analytic  results.  For  more  general  geometries  including  classically  chaotic  systems,  the  properties  of  the 
calculated  eigenvalue  spectrum,  such  as  the  energy  level-spacing  statistics,  agree  well  with  the  known  results 
for  different  symmetry  classes.  In  fact,  our  method  is  capable  of  finding  eigenstates  of  Dirac  fermions  under 
arbitrarily  electrical  potential  profiles.  Our  mafrix  formulation  can  be  applied  directly  to  one-dimensional 
systems.  Through  a  straightforward  extension  of  the  Dirac  spinor  to  four  components  and  a  proper  revision 
in  the  discretization  and  boundary  constraints,  the  method  can  be  extended  to  solving  the  Dirac  equation  in 
three  dimensions  as  well.  Our  main  finding  was  fhaf  the  relativistic,  spinor  type  of  wavefunctions  associated 
with  Dirac  fermions  can  be  highly  non-uniform  in  chaotic  billiards,  and  truly  relativistic  quantum  scars  do 
exist.  Examples  of  such  scars  are  shown  in  Figs.  1  and  2. 

A  number  of  areas  in  physics  can  benefit  from  an  efficient  method  for  solving  the  Dirac  equation.  The 
most  relevant  area  is  graphene  physics.  Graphene  ribbons  exhibit  a  linear  energy-momentum  relation  near  a 
Dirac  point  in  the  energy-band  diagram,  which  is  a  characteristic  of  relativistic  quantum  motion  of  massless 
fermions.  In  the  presence  of  short-range  potentials,  two  Dirac  points  are  coupled  together.  It  is  thus  of  basic 
interest  to  investigate  the  behavior  of  pure  Dirac  fermions  to  distinguish  them  from  those  due  to  the  coupling 
of  two  relativistic  particles. 

Details  of  this  work  can  be  found  in 

•  X.  Ni,  L.  Huang,  Y.-C.  Lai,  and  C.  Grebogi,  “Scarring  of  Dirac  fermions  in  chaotic  billiards,”  Physical 
Review  E  86,  015702,  1-8  (2012). 


3.1.2  Chiral  scars  in  chaotic  Dirac  fermion  systems 

During  the  project  period,  we  discovered  a  new  class  of  quantum  scars:  chiral  scars  in  chaotic  Dirac  fermion 
systems.  This  was  motivated  by  the  following  question:  are  there  characteristics  of  relativistic  quantum  scars 
which  differ  fundamentally  from  those  associated  with  non-relativistic  quantum  scars?  To  make  possible  our 
search  for  such  characteristics,  we  developed  an  analytic  approach  to  calculating  the  eigenstates  of  massless 
Dirac  fermions  in  a  broad  class  of  chaotic  billiards  by  using  the  method  of  conformal  mapping.  In  particular, 
for  any  shape  in  the  class,  a  proper  conformal  mapping  can  transform  it  to  a  shape  for  which  the  solutions 
of  the  Dirac  equation  can  be  written  down  analytically.  An  inverse  transform  of  the  solutions  thus  leads  to 
eigenstates  in  the  original  domain.  This  method  allows  us  to  calculate  an  unprecedentedly  large  number  of 
eigenvalues  and  eigenstates  with  high  accuracy.  Taking  advantage  of  this  powerful  method,  we  succeeded 
in  identifying  one  such  characteristic  associated  with  the  phase  of  the  wavefunction.  To  explain  this,  we 
note  that,  in  non-relativistic  quantum  systems,  when  a  particle  traverses  one  cycle  along  a  scarred  orbit, 
the  associated  quantum  phase  change  is  zero  or  2;i.  However,  when  the  various  eigenstate  solutions  of  the 
massless  Dirac  equation  were  examined,  we  found  one  subclass  of  scarred  orbits  for  which  one  complete 
itinerary  brings  about  a  phase  change  of  only  Jt.  In  fact,  it  takes  two  cycles  for  the  phase  of  the  wavefunction 
to  become  2jt  so  that  the  wavefunction  to  return  completely  to  its  original  value.  This  relativistic  quantum 
phenomenon  is  originated  from  the  chirality  of  the  massless  Dirac  fermions,  and  consequently  we  named 
such  scars  chiral  scars.  It  should  be  noted  that,  despite  the  emergence  of  chiral  scars,  majority  of  the 
scars  are  conventional  in  the  sense  that  the  phase  change  associated  with  one  cycle  is  2%.  We  developed  a 
semiclassical  theory  to  explain  the  physical  origin  of  chiral  scars,  which  are  uniquely  relativistic  quantum 
scars  and  find  no  counterparts  in  non-relativistic  quantum  systems. 

Details  of  this  work  can  be  found  in 

•  H.-Y.  Xu,  L.  Huang,  Y.-C.  Lai,  and  C.  Grebogi,  “Chiral  scars  in  chaotic  Dirac  fermion  systems,” 
Physical  Review  Letters  110,  064102,  1-5  (2013).  This  work  was  selected  by  PRL  Editors  as  an 
Editors’  Suggestion. 
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Figure  3:  Tunneling  rates  and  LDS  pat¬ 
terns  for  massless  Dirac  fermion  in  inte- 
grable  and  chaotic  double-well  systems, 

where  (|)  and  ^  are  two  components  of  the 
Dirac  spinor.  A  unit  system  h  =  c  =  \  has 
been  used  in  our  calculations,  in  which  the 
barrier  height  Vb  is  about  60.  For  the  rect¬ 
angle  double-well,  the  width  and  the  height 
are  set  to  2  and  1,  respectively.  The  bow- 
tie  chaotic  shape  is  obtained  by  cutting  the 
rectangle  with  three  arcs  such  that  the  cut 
parts  are  all  0.3  measured  on  both  sides  of 
the  baseline  and  central  vertical  line  of  the 
rectangle.  The  theoretical  ratio  of  the  left 
well  width  to  the  barrier  width  is  24: 1. 

^ '  '0 

3.1.3  Relativistic  quantum  tunneling  in  Dirac  fermion  and  graphene  systems 

In  2011,  the  remarkable  phenomenon  of  chaos-regularized  tunneling  in  non-relativistic  quantum  systems 
was  uncovered,  where  classical  chaos  can  suppress,  significantly,  the  spread  in  the  tunneling  rate  commonly 
seen  in  systems  whose  classical  dynamics  are  regular.  For  example,  for  a  system  consisting  of  two  symmet¬ 
rical  cavities  connected  by  a  one-dimensional  potential  barrier  along  the  line  of  symmetry,  when  the  classical 
dynamics  in  each  cavity  is  integrable,  for  sufficiently  large  energy  the  tunneling  rate  can  assume  many  val¬ 
ues  in  a  wide  interval.  Choosing  the  geometry  of  the  cavity  such  that  the  classical  dynamics  become  chaotic 
can  greatly  enhance  and  regularize  quantum  tunneling.  Heuristically,  this  can  be  understood,  as  follows. 
When  the  potential  barrier  is  infinite,  each  cavity  is  a  closed  system  with  an  infinite  set  of  eigenenergies 
and  eigenstates.  Many  eigenstates  are  concentrated  on  classical  periodic  orbits,  forming  quantum  scars.  For 
classically  integrable  cavity,  some  stable  or  marginally  stable  periodic  orbits  can  persist  when  the  potential 
barrier  becomes  finite  so  that  each  cavity  system  is  effectively  an  open  quantum  system.  Many  surviving 
eigenstates  correspond  to  classical  periodic  orbits  whose  trajectories  do  not  encounter  the  potential  barrier, 
generating  extremely  low  tunneling  rate  even  when  the  energy  is  comparable  with  or  larger  than  the  height 
of  the  potential  barrier.  The  eigenstates  corresponding  to  classical  orbits  that  interact  with  the  potential 
barrier,  however,  can  lead  to  relatively  strong  tunneling.  In  a  small  energy  interval  the  quantum  tunneling 
rate  can  thus  spared  over  a  wide  range.  Flowever,  when  the  classical  dynamics  is  chaotic,  isolated  orbits  that 
do  not  interact  with  the  potential  barrier  are  far  less  likely  and,  consequently,  the  states  associated  with  low 
tunneling  rates  disappear,  effectively  suppressing  the  spread  in  the  tunneling  rate. 

We  addressed  the  question  of  whether  chaos  can  regularize  tunneling  in  relativistic  quantum  systems.  To 
be  concrete,  we  studied  the  motion  of  massless  Dirac  fermions  in  the  setting  of  resonant  tunneling  to  facili¬ 
tate  comparison  with  the  non-relativistic  quantum  case.  Our  extensive  computations  revealed  unequivocally 
the  existence  of  surviving  eigenstates  that  lead  to  extremely  low  tunneling  rates.  As  for  the  non-relativistic 
quantum  case,  making  the  cavities  classically  chaotic  can  greatly  regularize  the  quantum  tunneling  dynam¬ 
ics.  Figure  3  shows  the  generalized  tunneling  rate  R  versus  the  normalized  energy  E /Vq  for  a  massless  Dirac 
fermion  in  the  double-well  barrier  system  for  two  types  of  geometry:  one  classically  integrable  and  another 
chaotic.  For  the  integrable  geometry,  we  observe  the  existence  of  states  with  extremely  low  tunneling  rates, 
as  indicated  by  the  arrow  in  Fig.  3(a).  These  correspond  to  states  localized  nearly  entirely  in  the  left  or  right 
side  of  the  potential  barrier,  which  “survive”  the  tunneling  process  between  the  two  sides,  as  indicated  by 
the  accompanying  pattern  of  local  density  of  states  (LDS).  In  non-relativistic  quantum  transport,  these  are 
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Figure  4:  liinneling  rates  and  LDS 
patterns  in  an  integrable  and  a  chaotic 
graphene  douhle-well  systems,  where  A 
and  B  denote  the  two  distinct  types  of 
atoms  in  a  graphene  unit  cell.  Tight- 
binding  Hamiltonian  is  used.  The  barrier 
height  is  fixed  at  Vb  ~  0.67t  (with  t  denot¬ 
ing  the  nearest-neighbor  hopping  integral). 
After  converting  to  the  /i  =  c  =  1  unit  sys¬ 
tem,  the  geometric  measures  and  the  bar¬ 
rier  potential  are  the  same  as  for  the  single 
Dirac  fermion  system. 


effectively  quantum  pointer  states.  In  relativistic  quantum  systems,  we  observe  that  both  components  of 
the  underlying  Dirac  spinor  exhibit  a  heavy  probability  concentration  on  orbits  along  which  particles  travel 
vertically  back  and  forth  on  either  side,  parallel  to  the  barrier.  For  the  chaotic  geometry,  while  signatures  of 
pointer  states  can  still  be  found,  they  are  much  weaker  than  those  in  the  integrable  counterpart,  as  shown  in 
Fig.  3(b)  and  the  accompanying  LDS  pattern. 

For  practical  implications,  we  considered  resonant  tunneling  devices  made  entirely  of  graphene  and 
calculated  the  tunneling  rate  for  different  energy  values.  We  obtained  qualitatively  similar  results  to  those 
for  massless  Dirac  fermions,  as  shown  in  Fig.  4.  One  unique  feature  for  both  the  Dirac  and  graphene  systems, 
which  finds  no  counterpart  in  non-relativistic  quantum  tunneling  devices,  is  the  high  tunneling  rate  in  the 
regime  where  the  particle  energy  is  smaller  than  the  height  of  the  potential  barrier.  This  is  a  manifestation 
of  Klein  tunneling.  To  explain  the  numerical  findings,  we  developed  a  theory  based  on  the  concept  of  self 
energies  and  the  complex  energy  spectrum  of  the  non-Hermitian  Hamiltonian  for  the  “open”  cavity. 

Details  of  this  work  can  be  found  in 

•  X.  Ni,  L.  Huang,  Y.-C.  Lai,  and  L.  M.  Pecora,  “Effect  of  chaos  on  relativistic  quantum  tunneling,” 
Europhysics  Letters  98,  50007,  1-6  (2012).  This  paper  was  highlighted  in  Europhysics  News,  Vol.  43, 
No.  6  (2012). 


3.1.4  Relativistic  quantum  tunneling  in  nonhyperbolic  chaotic  systems 

A  challenging  and  fundamental  problem  in  the  field  of  quantum  chaos  is  to  understand  quantum  manifesta¬ 
tions  of  chaotic  behaviors  in  situations  where  the  classical  dynamics  are  nonhyperbolic.  In  classical  Hamil¬ 
tonian  systems,  nonhyperbolicity  is  characterized  by  coexistence  of  chaos  and  Kolmogorov-Arnold-Moser 
(KAM)  tori  in  the  phase  space,  hence  the  term  “mixed  phase  space.”  The  problem  had  been  addressed  previ¬ 
ously  but  mostly  in  the  realm  of  nonrelativistic  quantum  mechanics  governed  by  the  Schrodinger  equation. 
For  example,  in  open  Hamiltonian  systems  where  the  classical  dynamics  are  chaotic  scattering,  presence  of 
KAM  tori  can  result  in  enhanced  fluctuations  of  the  semiclassical  S-matrix  elements  or  even  lead  to  fractal 
fluctuations. 

Development  of  graphene  physics  in  the  past  decade  stimulated  a  great  deal  of  interest  in  relativistic 
quantum  mechanics  as  applied  to  solid-state  devices.  Graphene  systems,  such  as  quantum  dots  that  exhibit 
chaos  in  the  classical  limit,  were  studied  with  respect  to  issues  such  as  the  energy-level  statistics,  quantum 
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scars,  and  scattering  dynamics.  From  the  classical  point  of  view,  nonhyperbolic  systems  are  generic  while 
integrable  and  fully  chaotic  systems  are  the  exception,  as  the  latter  correspond  to  the  two  opposite  extreme 
cases  in  the  spectrum  of  nonlinear  Hamiltonian  systems.  Thus,  in  applications  of  graphene  devices,  nonhy¬ 
perbolic  dynamics  are  expected  to  be  more  common  than  integrable  or  chaotie  dynamics,  and  this  demands 
a  good  understanding  of  the  relativistie  quantum  manifestations  of  nonhyperbolicity. 

The  linear  energy-momentum  relation  in  graphene  systems,  which  is  characteristic  of  relativistic  quan¬ 
tum  motion,  however,  holds  only  approximately  near  a  Dirac  point  in  the  energy-band  diagram.  In  addition, 
while  the  motion  of  electrons  in  graphene  near  a  single  Dirac  point  can  be  described  by  the  Dirac  equation, 
in  realistic  systems  the  coupling  between  motions  associated  with  the  coexisting  Dirac  points  cannot  be  ne¬ 
glected.  Thus,  strictly  speaking,  study  of  chaotic  graphene  systems  can  lead  only  to  a  partial  understanding 
of  various  phenomena  in  the  emerging  field  of  relativistic  quantum  chaos.  To  obtain  a  complete  picture, 
fermion  systems  governed  by  the  Dirac  equation  must  be  studied. 

A  paradigm  of  nonhyperbolic  chaotic  systems  is  the  “mushroom”  billiard  for  which  a  mathematical 
proof  of  nonhyperbolicity  was  obtained  by  Bunimovich  in  2001.  By  placing  a  hnite  potential  barrier  along 
the  vertical,  symmetric  line  of  the  billiard,  we  effectively  generate  the  setting  of  resonant  tunneling  where  a 
particle  can  tunnel  through  the  barrier  from  one  side  of  the  billiard  to  the  other.  We  assumed  that  the  par¬ 
ticle  is  a  massless  Dirac  fermion,  whose  motion  is  governed  by  the  Dirac  equation,  and  we  focused  on  the 
effect  of  nonhyperbolicity  on  the  tunneling  rate  y{E),  the  fundamental  quantity  characterizing  the  quantum 
tunneling  dynamics,  which  depends  on  the  particle  energy  E.  Our  main  finding  was  that  nonhyperbolicity 
leads  to  a  “clustering”  phenomenon  where  the  majority  of  the  values  of  y{E)  fall  into  a  narrow  “band”  in 
the  {E,y)  plane,  with  relatively  few  values  outside  the  band,  as  shown  in  Fig.  5.  To  understand  this  cluster¬ 
ing  phenomenon,  we  develop  a  theoretical  framework  based  on  the  concept  of  self-energies  to  calculate  the 
tunneling  rate.  We  then  analytically  solved  the  Dirac  equation  both  in  one  dimension  and  in  two  dimensions 
for  a  circular-ring  type  of  tunneling  system  with  integrable  dynamics  in  the  classical  limit.  Due  to  the  pres¬ 
ence  of  relatively  few  and  distinct  classical  periodic  orbits  in  the  integrable  component,  the  corresponding 
relativistic  quantum  states  can  have  drastically  different  behaviors,  leading  to  a  wide  spread  of  the  values  of 
the  tunneling  rate  in  the  {E,y)  plane.  However,  the  chaotic  component  provides  significantly  more  quantum 
states  for  tunneling  due  to  the  inhnite  set  of  unstable  periodic  orbits  embedded  therein.  These  states  are 
dynamically  “similar”  because  of  the  ergodic  nature  of  chaos,  and  this  leads  to  the  clustering  of  the  values 
of  the  tunneling  rate  in  a  narrow  band. 

Details  of  this  work  can  be  found  in 

•  X.  Ni,  L.  Huang,  L.  Ying,  and  Y.-C.  Lai,  “Relativistic  quantum  tunneling  of  Dirac  fermions  in  nonhy¬ 
perbolic  chaotic  systems,”  Physical  Review  B  87,  224304,  1-10  (2013). 


3.1.5  Effect  of  many-body  interactions  on  relativistic  quantum  chaos 

In  quantum  chaos,  the  most  often  studied  setting  was  that  of  single-particle  quantum  dynamics.  While 
there  were  previous  studies  of  the  interplay  between  many-body  interactions  and  classical  chaos,  these  were 
exclusively  for  non-relativistic  quantum  systems  described  by  the  Schrddinger  equation.  To  investigate  the 
effect  of  chaos  on  relativistic  quantum  systems  with  many-body  interactions  is  an  outstanding  problem,  yet 
it  is  not  only  fundamental  to  physics,  but  also  important  for  the  practical  development  of  relativistic  quantum 
devices  based  on  graphene  or  alternative  two-dimensional  Dirac  material  systems. 

In  the  project  period,  we  studied  quantum  chaos  in  the  presence  of  many-body  interactions  by  using  the 
standard  Hubbard  model.  This  paradigmatic  model  to  treat  interacting  particles  in  a  lattice  was  originally 
proposed  in  1963  to  describe  the  transition  between  conducting  and  insulating  systems.  For  electrons  in  a 
solid,  comparing  with  the  conventional  tight-binding  model  representing  a  single  electron  Hamiltonian,  the 
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Figure  5:  Tunneling  rates  versus  the  en¬ 
ergy  and  representative  eigenstates  for 
mushroom  billiard  made  of  graphene. 

The  dimensions  of  the  graphene  system 
and  the  barrier  height  and  width  are  pro¬ 
portional  to  those  of  the  nanoscale  Dirac- 
fermion  tunneling  system  that  we  studied, 
but  the  first  and  the  second  rows  of  pat¬ 
terns  represent,  instead  of  Dirac  spinors, 
the  A  and  B  atoms  in  the  graphene  unit  cell. 
The  eigenenergies  and  eigenstates  are  cal¬ 
culated  by  using  the  standard  tight-binding 
Hamiltonian  for  the  corresponding  closed 
graphene  system. 


Figure  6:  Effect  of  many-body  in¬ 
teractions  on  eigenstates.  Representa¬ 
tive  confined  eigenstates  associated  with 
spin-up  (red/solid  curves)  and  spin-down 
(blue/dashed  curves)  electrons:  (a)  prob¬ 
ability  density  profile  of  one  eigenstate 
without  a  potential  barrier,  where  it  extends 
in  both  potential  wells,  (b)  probability  den¬ 
sity  profile  of  the  eigenstate  with  a  poten¬ 
tial  barrier  at  .r  =  0  (the  gray  rectangle).  In 
(b),  there  is  spin  polarization,  i.e.,  spin-up 
electrons  reside  in  the  right  well  and  spin- 
down  electrons  reside  in  the  left  well.  The 
insets  in  both  panels  show  the  correspond¬ 
ing  patterns  of  local  density  of  state  in  the 
entire  domain. 


Hubbard  model  contains  a  potential  term  to  include  the  many-body  effect  through  the  mechanism  of  on-site 
Coulomb  interaction.  There  has  been  a  great  deal  of  interest  in  the  Hubbard  model  due  to  its  relevance  to 
frontier  problems  in  condensed  matter  physics  such  as  high-temperature  superconductivity  and  the  trapping 
of  ultra-cold  atoms  in  optical  lattices.  While  the  Hubbard  model  is  much  more  challenging  and  sophisticated 
than  the  tight-binding  model,  it  can  serve  as  a  paradigm  to  gain  signihcant  physical  insights  into  many-body 
relativistic  quantum  manifestations  of  distinct  type  of  classical  dynamics. 

In  our  research,  we  focused  on  graphene  systems  and  studied  the  particular  phenomenon  of  quantum 
resonant  tunneling.  The  typical  setting  of  a  quantum  tunneling  system  consists  of  two  symmetric  potential 
wells  separated  by  a  potential  barrier  in  between.  The  whole  system,  which  includes  the  left  and  right  wells 
as  well  as  the  barrier,  is  closed,  and  its  geometrical  shape  can  he  chosen  to  yield  characteristically  distinct 
types  of  dynamics  in  the  classical  limit. 

In  order  to  uncover  the  unique  relativistic  quantum  phenomena  caused  hy  classical  chaos  in  the  presence 
of  many-body  interactions,  we  hrst  studied  the  class  of  integrable  systems  of  rectangular  shape.  Since 
the  whole  system  is  closed,  we  calculated  the  eigenenergies  and  investigated  various  eigenstates  from  the 
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mean-field  Hubbard  Hamiltonian.  A  striking  finding  is  the  emergence  of  a  class  of  eigenstates  with  near 
zero  tunneling.  In  particular,  for  such  an  eigenstate,  the  spin-up  and  spin-down  wavefunctions  are  typically 
separated,  i.e.,  the  spin-up  electrons  reside  in  only  one  potential  well  while  the  spin-down  electrons  reside 
in  the  other,  as  exemplified  in  Fig.  6.  As  a  result,  if  the  initial  state  is  spin-up  in  one  potential  well,  it 
is  localized  and  will  stay  in  the  same  well  practically  for  an  inhnite  amount  of  time  with  little  quantum 
tunneling.  When  the  potential  term  describing  the  on-site  Coulomb  interactions  is  removed  so  that  the 
Hamiltonian  becomes  that  of  the  tight-binding  type,  such  localized  states  no  longer  exist,  indicating  that 
they  are  result  of  electron-electron  interactions  and  consequently  a  distinct  many-body  phenomenon.  We 
derived  an  approximate  theory,  based  on  the  simplified  picture  of  one-dimensional  tunneling  of  massless 
Dirac  fermions,  to  explain  the  physical  origin  of  the  localized  states.  We  further  found  that,  when  the 
geometrical  shape  is  that  of  stadium  or  bowtie  so  that  the  classical  dynamics  is  chaotic,  the  localized  states 
are  greatly  suppressed  and  the  tunneling  rates  are  signihcant.  This  means  that,  classical  chaos  is  capable  of 
effectively  removing  the  localized  states.  In  addition  to  the  classically  integrable  and  fully  chaotic  domains, 
we  also  considered  a  class  of  domains,  the  mushroom  billiard,  in  which  the  classical  dynamics  is  mixed  or 
nonhyperbolic  with  coexisting  regular  and  chaotic  components  in  the  phase  space.  We  found  that,  due  to  the 
chaotic  component,  quantum  tunneling  can  also  be  regularized  and  enhanced. 

From  the  standpoint  of  device  application,  the  localized  states  present  an  obstacle  to  effective  tunneling 
and  such  states  are  therefore  undesirable.  The  main  message  is  then  that  chaos  can  signihcantly  enhance 
the  tunneling  process  in  realistic  situations  where  electron-electron  interactions  are  present.  This  implies 
that  classical  chaos  is  capable  of  facilitating  greatly  relativistic  quantum  tunneling,  which  is  desirable  in  the 
development  of  nanoscale  devices  such  as  graphene-based  resonant-tunneling  diodes. 

Details  of  this  work  can  be  found  in 

•  L.  Ying,  G.-L.  Wang,  L.  Huang,  and  Y.-C.  Lai,  “Quantum  chaotic  tunneling  in  graphene  systems  with 
electron-electron  interactions,”  Physical  Review  B  90,  224301,  1-14  (2014). 

3.2  Conductance  fluctuations  associated  with  quantum  transport  in  graphene  systems 

A  fundamental  problem  in  quantum  transport  through  nanoscale  devices  is  conductance  fluctuations.  Con¬ 
sider,  for  example,  a  quantum-dot  system.  As  the  Fermi  energy  of  the  conducting  electrons  is  varied,  the 
conductance  can  exhibit  fluctuations  of  distinct  characteristics,  depending  on  the  geometrical  shape  of  the 
dot.  Research  in  the  past  two  decades  demonstrated  that  the  nature  of  the  corresponding  classical  dynamics 
can  play  a  key  role  in  the  conductance-fluctuation  pattern.  For  example,  when  the  classical  scattering  dy¬ 
namics  is  integrable  or  has  a  mixed  phase-space  structure,  there  can  be  sharp  fluctuations  in  the  conductance 
curve.  However,  when  the  classical  dynamics  is  fully  chaotic,  the  conductance  variations  tend  to  be  more 
smooth. 

Qualitatively,  why  the  conductance-fluctuation  patterns  depend  on  the  nature  of  the  corresponding  clas¬ 
sical  dynamics  can  be  understood,  as  follows.  Sharp  conductance  fluctuations  are  caused  by  quantum  pointer 
states,  which  are  resonant  states  of  hnite  but  long  lifetime  formed  inside  the  nanostructure,  a  kind  of  Fano 
resonance.  For  example,  for  a  quantum-dot  system  whose  classical  dynamics  is  regular  or  contains  a  sig¬ 
nificantly  regular  component,  there  are  stable  periodic  orbits.  For  the  corresponding  closed  dot  geometry, 
highly  localized  states  can  form  about  the  periodic  orbits,  forming  quantum  scars.  When  leads  are  attached 
to  the  quantum  dot  so  that  the  system  becomes  open,  some  periodic  orbits  can  still  survive,  leading  to  res¬ 
onant  states,  or  quantum  pointer  states.  Since  the  corresponding  classical  orbits  are  stable,  the  resonant 
states  can  have  long  lifetime,  so  they  couple  to  the  leads  only  weakly.  As  a  result,  very  narrow  resonances 
can  form  about  the  energy  values  that  are  effectively  the  eigenenergies  in  the  corresponding  closed  system. 
When  a  modification  to  the  dot  geometry  is  introduced  so  that  the  underlying  classical  dynamics  becomes 
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fully  chaotic,  no  stable  periodic  orbits  can  exist.  While  scars  can  still  be  formed  around  classically  unstable 
periodic  orbits  in  the  closed  system,  the  corresponding  resonant  states  in  the  open  system  generally  will 
have  much  shorter  lifetimes,  effectively  eliminating  the  narrow  resonances  in  the  conductance  fluctuation 
pattern.  This  understanding  also  led  to  the  idea  that  the  interplay  between  the  classical  dynamics  and  the 
conductance-fluctuation  pattern  can  be  exploited  for  effectively  modulating  quantum  transport  (Sec.  3.2.2). 

3.2.1  Insensitivity  of  conductance  fluctuations  to  the  nature  of  classical  dynamics 

We  studied  conductance  fluctuation  patterns  in  graphene  devices.  In  a  graphene  quantum  dot,  the  character¬ 
istics  of  conductance  fluctuations  depend  on  the  nature  of  the  classical  dynamics  in  a  way  similar  to  that  for 
conventional  two-dimensional  semiconductor  quantum  dot  systems.  However,  the  magnetic  properties  of 
graphene  are  quite  different  from  those  associated  with  conventional  two-dimensional  electrical  gas  (2DEG) 
systems  because  of  the  special  lattice  structure  of  graphene  and  the  relativistic  quantum  behavior  of  quasi¬ 
particles.  For  example,  in  graphene  quantum  Hall  effect  can  be  observed  even  at  room  temperature  due  to 
the  massless  Dirac  fermion  nature  of  quasiparticles  and  significantly  reduced  scattering  effects.  Especially, 
the  linear  energy-momentum  relation  in  graphene  stipulates  that  the  Eandau  levels  are  distributed  according 
to  ±\/iV,  where  N  is  the  Landau  index,  versus  the  proportional  dependence  on  N  in  conventional  2DEG 
systems. 

Our  primary  goal  was  to  uncover  any  unique  characteristics  of  conductance  fluctuations  in  presence  of 
a  magnetic  field,  in  fhe  framework  of  non-equilibrium  Green’s  funclion  (NEGF).  Our  investigafion  revealed 
fwo  surprising  phenomena.  Firsf,  in  the  parameter  plane  spanned  by  the  strength  of  a  perpendicular  mag¬ 
netic  field  and  fhe  Fermi  energy,  there  are  regions  of  regular  and  random  conductance  oscillations.  As  the 
Fermi  energy  or  the  magnetic  field  strengfh  is  changed,  the  fluctuations  can  be  quite  regular  and  then  ran¬ 
dom,  and  vice  versa,  implying  a  kind  of  “coexistence”  of  regular  and  irregular  conductance  fluctuations  as 
a  single  physical  parameter  is  varied.  Second,  there  exists  a  universal  scaling  behavior  characterizing  the 
magnetic  conductance-fluctuation  patterns,  which  does  not  depend  on  the  nature  of  the  classical  dynamics. 
Through  examining  the  edge  states  between  adjacent  Landau  levels  and  the  correspondence  of  their  for¬ 
mation  to  conductance  resonances,  we  found  that,  as  the  energy  is  increased  or  the  magnetic-field  sfrength 
is  decreased  so  thaf  a  Landau  level  is  crossed,  a  new  set  of  fluctuation  patterns  emerges,  leading  to  nearly 
periodic  oscillations  but  of  different  period  in  the  conductance  versus  the  Fermi  energy  or  the  magnetic-field 
strengfh.  Irregularity  associated  with  the  conductance  fluctuations  is  developed  as  more  periodic  patterns  set 
in.  Considering  that  previous  works  pointed  to  the  critical  role  played  by  the  nature  of  classical  dynamics  in 
conductance  fluctuations,  our  finding  fhaf  fhe  presence  of  magnefic  field  can  greatly  suppress  this  sensitivity 
to  classical  dynamics  was  striking. 

One  pertinent  question,  when  a  perpendicular  magnetic  field  is  applied  to  a  quantum  dot,  is  how  conduc¬ 
tance  fluctuations  are  affected  by  the  size  of  the  quantum  dot.  In  a  previous  experimental  study  of  quantum 
dots  of  size  ranging  from  0.7  /urn  to  1.2  qm,  nearly  periodic  conductance  oscillations  had  been  found  as 
the  magnetic-field  strength  is  varied.  The  frequency  of  the  oscillation  pattern,  the  so-called  magnetic  fre¬ 
quency,  follows  a  scaling  relation  with  the  edge  size  of  the  dot.  In  a  study  of  the  magnetic  scaling  behavior 
for  graphene  quantum  dots,  it  was  found  that  for  small  dots  of  edge  size  less  than  0.3  qm,  the  magnetic 
frequency  exhibits  a  scaling  relation  with  the  dot  area.  We  focused  on  a  particular  set  of  scarred  orbits 
and  examined  the  resulting  conductance  oscillations.  We  found,  for  graphene  quantum  dots,  below  the  first 
Landau  level,  the  conductance  exhibits  periodic  oscillations  with  the  magnetic-field  sfrength  and  with  the 
Fermi  energy.  The  magnetic  frequency  scales  linearly  with  the  dot  size.  However,  the  energy  frequency,  the 
inverse  of  the  variation  in  the  Fermi  energy  for  the  conductance  to  complete  one  cycle  of  oscillation,  scales 
inversely  with  the  dot  size.  Beyond  the  regime  of  periodic  conductance  oscillations,  new  sets  of  scarred 
orbits  emerge  and  evolve,  each  with  its  own  period,  leading  to  random-like  conductance  fluctuations.  The 
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Figure  7:  A  possible  experimental 
scheme  to  harness  transport  through 
a  semiconductor  2DEG  quantum-dot 

system,  where  2DEG  is  formed  at  the 
Q  GaAs/Alo.sGaojAs  hetero-interface.  The 
heterostructure  sits  on  a  n+Si  substrate 
(purple),  covered  by  300  nm  Si02  (blue) 
and  contacted  by  Au/Cr  (yellow).  When 
a  suitable  gate  voltage  is  applied  to  gener¬ 
ate  a  circular  forbidden  region  (for  classi¬ 
cal  orbits)  at  the  center  of  the  device,  the 
resulting  closed  system  is  a  Sinai  billiard. 

Open  quantum-dot  system  can  be  formed 
by  attaching  leads  to  the  billiard  system.  In 
our  study,  we  placed  the  leads  in  the  mid¬ 
dle  of  the  dot,  as  shown  in  the  second  row. 
Similar  idea  can  be  applied  to  graphene  to 
implement  the  chaos-based  harnessing  of 
quantum  transport. 

remarkable  feature  is  that  these  scaling  behaviors  are  independent  of  the  nature  of  the  underlying  classical 
dynamics,  i.e.,  regular  or  chaotic. 

Details  of  this  work  can  be  found  in 

•  L.  Ying,  L.  Huang,  Y.-C.  Lai,  and  C.  Grebogi,  “Conductance  fluctuations  in  graphene  systems:  The 
relevance  of  classical  dynamics,”  Physical  Review  B  85,  245448,  1-7  (2012). 


3.2.2  Harnessing  quantum  transport  by  transient  chaos 

Controlling  chaos  in  dynamical  systems  has  been  studied  for  more  than  two  decades  since  the  seminal  work 
of  Ott,  Grebogi,  and  Yorke  in  1990.  The  basic  idea  was  that  chaos,  while  signifying  random  or  irregular 
behavior,  should  not  be  viewed  as  a  nuisance  in  applications  of  nonlinear  dynamical  systems.  In  fact,  given 
a  chaotic  system,  that  there  are  an  infinite  number  of  unstable  periodic  orbits  embedded  in  the  underlying 
chaotic  invariant  set  means  that  there  are  an  equally  infinite  number  of  choices  for  the  operational  state 
of  the  system  depending  on  need,  provided  that  any  such  state  can  be  stabilized.  Then,  the  intrinsically 
sensitive  dependence  on  initial  conditions,  the  hallmark  of  any  chaotic  system,  implies  that  it  is  possible  to 
apply  small  perturbations  to  stabilize  the  system  about  any  desirable  state.  Controlling  chaos  has  since  been 
studied  extensively  and  examples  of  successful  experimental  implementation  abound  in  physical,  chemical, 
biological,  and  engineering  systems.  The  vast  literature  on  controlling  chaos,  however,  was  limited  to 
nonlinear  dynamical  systems  in  the  classical  domain.  We  elaborated  the  idea  that  chaos  may  be  exploited  to 
harness  or  modulate  quantum-mechanical  behaviors  by  conducting  a  systematic  study  of  the  role  of  transient 
chaos  in  modulating  quantum  transport  dynamics.  In  particular,  in  applications  such  as  the  development  of 
electronic  circuits  and  nanoscale  sensors,  severe  conductance  fluctuations  are  undesirable  and  are  to  be 
eliminated  so  that  stable  device  operation  can  be  achieved.  An  outstanding  question  was  then,  can  practical 
and  experimentally  feasible  schemes  be  articulated  to  modulate  the  quantum  conductance  fluctuations?  We 
demonstrated  and  provided  theoretical  understanding  that  classical  transient  chaos  can  be  used  to  effectively 
modulate  conductance-fluctuation  patterns  associated  with  quantum  transport  through  nanostructures. 

A  quantum  dot  typically  consists  of  a  finite  device  region  of  certain  geometrical  shape,  such  as  a  square, 
a  circle,  or  a  stadium,  and  a  number  of  leads  connected  with  the  device  region.  To  demonstrate  quantum 
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Figure  8:  Exploiting  chaos  for  quan¬ 
tum  control.  (a)  Conductance  ver¬ 
sus  Fermi  energy  for  four  semiconduc¬ 
tor  2DEG  quantum-dot  systems  (bottom 
to  top):  rectangular  dot,  rectangular  dot 
with  a  rectangular  forbidden  region  of  area 
0.25  pmx  0.25  /jm,  and  Sinai  dots  of  radii 
R=0.14  /jin  and  R=0.28  /jm.  The  area  of 
the  original  rectangular  dot  is  1  pmx  1  /jin 
and  the  lead  width  is  0.22  pm.  (b)  Con¬ 
ductance  fluctuation  patterns  for  graphene 
quantum  dots  of  the  same  geometry  as 
in  (a).  The  curves  in  (a)  and  (b)  were 
shifted  vertically  by  some  arbitrary  values 
to  facilitate  visualization,  (c)  Normalized 
number  of  sharp  resonances,  which  is  de¬ 
fined  to  quantify  the  degree  of  transmis¬ 
sion  fluctuations,  corresponding  to  the  four 
cases  in  (a),  i.e.,  the  top  two  rectangular 
symbols  represent  two  rectangular  quan¬ 
tum  dots  while  the  bottom  two  circle  sym¬ 
bols  correspond  to  two  Sinai  quantum  dots. 


harnessing  by  using  chaos,  we  conceived  generating  a  region  about  the  center  of  the  device  or  structure 
with  high  potential  so  that  it  is  impenetrable  to  classical  particles.  For  example,  consider  a  rectangular 
quantum  dot.  When  it  is  closed,  the  corresponding  classical  dynamics  is  integrable  so  that  extremely  narrow 
resonances  can  arise  in  the  quantum  transport  dynamics  of  the  open-dot  system.  Now  imagine  applying 
a  gate  voltage  to  generate  a  circular,  classically  forbidden  region  about  the  center  of  the  dot,  as  shown 
schematically  in  Fig.  7.  In  general,  the  potential  profile  will  be  smooth  in  space.  However,  the  scattering 
behavior  is  qualitatively  similar  to  that  from  an  infinite  potential  well.  Thus  in  our  simulation  we  adopted 
the  infinite  potential-well  assumption  for  the  central  region,  which  defines  a  “forbidden”  region.  Varying 
the  voltage  Vq  can  change  the  effective  radius  R  of  the  forbidden  region.  Classically,  the  closed  system  is 
thus  a  Sinai  billiard,  which  is  fully  chaotic,  insofar  as  the  radius  of  the  central  potential  region  R  is  not  zero. 
When  leads  are  connected  to  the  device  region  so  as  to  open  the  system,  chaos  becomes  transient.  The 
dynamical  characteristics  of  the  underlying  chaotic  invariant  set  can  be  adjusted  in  a  continuous  manner  by 
increasing  the  radius  R.  Quantum  mechanically  we  thus  expect  to  observe  increasingly  smooth  variations  in 
the  conductance  with,  e.g.,  the  Fermi  energy,  which  we  demonstrated  using  both  semiconductor  2DEG  and 
graphene  systems,  as  shown  in  Fig.  8. 

Details  of  this  work  can  be  found  in 

•  R.  Yang,  L.  Huang,  Y.-C.  Lai,  C.  Grebogi,  and  L.  M.  Pecora,  “Harnessing  quantum  transport  by 
transient  chaos,”  Chaos  23,  013125,  1-9  (2013). 


3.2.3  Quantum  chaotic  scattering  in  graphene  systems  without  invariant  classical  dynamics 

In  open  Hamiltonian  systems,  a  fundamental  type  of  chaotic  behaviors  is  transient  chaos,  which  leads  to 
chaotic  scattering.  Quantum  chaotic  scattering,  referred  to  as  the  study  of  the  unique  quantum  behaviors 
caused  by  chaotic  scattering  in  the  corresponding  classical  system,  is  highly  relevant  to  a  number  of  fields 
in  physics,  such  as  atomic  physics,  condensed  matter  physics,  and  acoustics. 

The  fundamental  quantity  characterizing  a  quantum  scattering  system  is  the  scattering  matrix,  or  the 
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S-matrix,  whose  elements  are  the  transition  probabilities  between  quantum  states  of  the  system  before  and 
after  the  scattering.  The  formulation  of  the  S-matrix  in  terms  of  classical  quantities  had  been  of  great 
interest  in  chemical  physics  even  before  chaos  started  attracting  wide  attention.  The  seminal  contribution 
by  Miller  in  1975,  who  obtained  a  formula  for  S-matrix  elements  in  terms  of  purely  classical  quantities 
in  the  semiclassical  regime  for  reactive  scattering  systems,  became  the  fundamental  tool  in  the  study  of 
quantum  chaotic  scattering.  Given  a  system  that  exhibits  chaotic  scattering  in  the  classical  limit,  the  S- 
matrix  elements  in  the  semiclassical  regime  exhibit  random  fluctuations  as  some  physical  parameters  of  the 
system,  such  as  the  energy  of  the  scattering  particle  or  the  strength  of  some  externally  applied  magnetic 
field,  change  in  a  classically  small  but  quantum-mechanically  large  range.  Depending  on  whether  classical 
scattering  is  hyperbolic  (fully  chaotic)  or  nonhyperbolic  (mixed),  the  statistical  properties  of  the  fluctuations 
in  the  S-matrix  elements  can  be  quite  distinct.  In  particular,  the  fluctuation  patterns  can  be  characterized  by 
the  energy-correlation  function,  where  a  faster  decay  of  the  function  from  unity  as  the  energy  difference  is 
increased  points  to  more  severe  fluctuation  patterns.  In  their  seminal  work  in  1988,  Bliimel  and  Smilansky 
showed  that  the  energy-correlation  function  is  related  to  the  decay  law  of  trajectories  in  the  classical  phase 
space  as  a  Fourier-transform  pair.  Thus,  if  classical  trajectories  decay  faster  from  the  scattering  region,  the 
energy-correlation  function  decreases  more  slowly  from  unity,  and  vice  versa. 

A  tacit  assumption  employed  in  previous  works  on  quantum  chaotic  scattering  is  that,  as  a  physical 
parameter  changes,  the  classical  dynamics  is  invariant  so  that  the  decay  law  remains  unchanged,  rendering 
meaningful  its  Fourier  transform.  Such  situations  arise  in  physical  systems,  for  example,  transport  of  elec¬ 
trons  through  a  quantum  dot  in  the  absence  of  magnetic  field  where  the  corresponding  classical  dynamics  is 
essentially  that  of  a  open  billiard  and  so  does  not  depend  on  the  electron  energy.  However,  when  a  magnetic 
field  is  present,  the  characteristics  of  classical  dynamics  can  change  drastically  with  the  particle  energy.  In 
particular,  in  an  early  work,  Breymann,  Kovacs  and  Tel  studied  chaotic  scattering  of  charged  particles  in 
an  open  three-disk  billiard  subject  to  a  perpendicular  magnetic  field,  which  is  effectively  a  three-terminal 
quantum-dot  system.  Classically,  since  the  Lorentz  force  depends  on  both  the  particle  energy  and  the  mag¬ 
netic  field,  even  for  a  fixed  magnetic  field  strength  the  dynamics  will  depend  on  the  energy.  We  thus  faced  a 
challenge  that,  as  the  particle  energy  is  systematically  changed,  there  is  no  unique  classical  correspondence. 
Quantum  mechanically,  the  scattering  matrix  elements  will  exhibit  fluctuations  with  the  particle  energy,  but 
the  Fourier-transform  formula  relating  the  energy  autocorrelation  function  and  classical  particle  decay  law 
is  no  longer  applicable.  How  could  then  the  fluctuation  patterns  be  characterized  and  understood? 

To  investigate  quantum  chaotic  scattering  in  absence  of  invariant  classical  dynamics,  we  used  quantum- 
transport  systems  in  two-dimensions  in  the  presence  of  a  perpendicular  magnetic  field  as  a  paradigm.  In 
particular,  utilizing  multi-terminal  graphene  quantum  dots  as  a  prototypical  class  of  systems,  we  articulated 
that  quantum  chaotic  scattering  can  be  physically  characterized  by  the  complex  eigenvalues  of  the  non- 
Hermitian  Hamiltonian  of  the  underlying  system.  A  quantum-dot  system  can  in  general  be  decomposed 
into  two  parts:  a  closed  device  or  scattering  region  and  the  set  of  semi-infinite  electronic  waveguides  (or 
leads).  The  device  Hamiltonian  is  Hermitian  and  permits  a  discrete  set  of  eigenvalues,  but  the  leads  possess 
a  continuum  of  energy  spectrum,  mathematically  described  by  complex  self-energy  terms.  Adding  the  self¬ 
energies  to  the  device  Hamiltonian  leads  to  a  non-Hermitian  Hamiltonian  describing  the  entire  open  quantum 
system,  whose  eigenvalues  are  in  general  complex.  The  key  physics  to  quantum  chaotic  scattering  is  that 
the  inverse  of  the  imaginary  part  of  a  complex  eigenvalue  is  nothing  but  the  lifetime  of  the  corresponding 
eigenstate  in  the  scattering  region.  Thus,  when  an  energy  eigenvalue  possesses  an  extremely  small  imaginary 
part,  the  corresponding  eigenstate  has  a  long  lifetime,  which  is  effectively  a  strongly  localized  state.  The 
coupling  between  such  a  localized  state  and  the  leads,  or  the  quantum  “environment”  of  the  device,  must 
necessarily  be  weak.  This  leads  to  an  abrupt  change  in  the  transmission  or  conductance  over  an  extremely 
small  energy  scale,  the  so-called  Fano  resonance.  Overall,  in  order  to  explain  a  numerically  calculated  or 
experimentally  observed  conductance-fluctuation  pattern  in  quantum-dot  systems  in  presence  of  a  magnetic 


14 


Geometry 


(C) 


Figure  9:  Conductance  fluctuations 
and  characterization  in  a  four-terminal 
graphene  quantum-dot  system  under  a 
perpendicular  magnetic  field:  (a)  geom¬ 
etry  of  the  system,  where  the  scattering 
(device)  region  is  denoted  by  thick  black 
curves  and  the  region  in  which  the  perpen¬ 
dicular  magnetic  field  is  applied  is  repre¬ 
sented  by  the  blue  dashed  circle;  (b)  val¬ 
ues  of  the  eigenenergies  in  the  complex 
plane,  where  eigenvalues  with  extremely 
small  imaginary  part  begin  to  appear  after 
the  first  Landau  level  arises;  (c)  a  represen¬ 
tative  strongly  localized  state  with  weak 
coupling  to  the  leads,  corresponding  to  the 
filled  red  circle  in  (b);  (d)  the  correspond¬ 
ing  local  conductance-fluctuation  pattern, 
where  the  vertical  red  dashed  line  corre¬ 
sponds  to  the  real  energy  value  of  the  filled 
red  circle  in  (b). 


field,  one  can  calculate  the  complex  eigenenergies  of  the  corresponding  non-Hermitian  Flamiltonian.  Energy 
regions  in  which  the  eigenenergies  exhibit  small  imaginary  parts  are  the  regions  where  severe  conductance 
fluctuations  can  be  anticipated.  These  results  should  provide  insights  into  the  development  of  graphene- 
based  quantum  transport  systems. 

Some  representative  results  are  shown  in  Fig.  9. 

Details  of  this  work  can  be  found  in 

•  G.-L.  Wang,  L.  Ying,  Y.-C.  Lai,  and  C.  Grebogi,  “Quantum  chaotic  scattering  in  graphene  systems  in 
the  absence  of  invariant  classical  dynamics,”  Physical  Review  E  87,  052908,  1-9  (2013). 


3.2.4  Universal  transform  for  Fano  resonance 

A  fundamental  phenomenon  associated  with  quantum  or  wave  scattering  dynamics  is  Fano  resonance.  A 
typical  scattering  system  consists  of  incoming  channels,  a  scatterer  or  a  conductor,  and  outgoing  channels. 
The  scatterer,  when  isolated,  can  be  regarded  as  a  closed  system  with  a  discrete  spectrum  of  intrinsic  energy 
levels.  When  the  energy  of  the  incoming  particle  or  wave  matches  an  energy  level,  a  resonant  behavior  can 
arise  in  some  experimentally  measurable  quantities,  such  as  the  conductance  in  a  quantum-transport  system. 
The  resonance  profile  is  typically  asymmetric,  and  can  in  general  be  expressed  as  (s  +  q)^ / (s^  +  1),  where 
s  is  the  normalized  energy  deviation  from  the  center  of  the  resonance,  and  q  is  the  parameter  characterizing 
the  degree  of  asymmetry  of  the  resonance.  The  asymmetry  parameter  q  is  of  great  experimental  importance 
as  it  determines  how  experimental  data  can  be  fitted  by  the  Fano  profile. 

The  profile  was  first  derived  by  Fano  in  the  study  of  inelastic  scattering  of  electrons  off  the  helium 
atom  and  auto-ionization,  although  the  phenomenon  was  predicted  earlier  in  elastic  neutron  scattering.  Be¬ 
ing  a  general  wave  interference  phenomenon,  Fano  resonance  has  been  found  in  many  contexts  in  physics, 
such  as  photon-ionization,  Raman  scattering,  photon-absorption  in  quantum-well  structure,  scanning  mi¬ 
croscopy  tunneling  in  the  presence  of  impurity,  transport  through  a  single-electron  transistor  and  through 
the  Ahamov-Bohm  interferometer,  transport  through  crossed  carbon  nanotubes,  microwave  scattering,  plas- 
monic  nanostructures  and  metamaterials,  and  optical  resonances.  Applications  exploiting  Fano  resonance 
have  even  been  proposed  for  biochemical  sensors. 
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An  unsettled  and  somewhat  controversial  issue  was  whether  the  asymmetric  parameter  q  should  be 
real  or  complex.  The  issue  is  of  high  experimental  relevance.  For  example,  it  was  demonstrated  that,  in 
single-channel  scattering  quantum-transport  systems,  q  is  strictly  real  if  the  time-reversal  symmetry  (TRS) 
is  preserved.  When  the  TRS  is  broken,  q  will  generally  be  complex,  which  was  observed  experimentally 
in  an  Aharnov-Bohm  interferometer  with  a  quantum  dot  embedded  in  one  of  the  arms.  With  a  changing 
magnetic  field,  q  oscillates,  which  was  explained  but  again  by  using  the  single-channel  scattering  model. 
Later,  it  was  shown  that  for  multi-channel  scattering,  q  is  in  general  complex  even  when  TRS  is  not  broken. 
It  was  also  proposed  that  complex  q  parameter  can  be  used  to  characterize  the  dephasing  time.  The  variation 
in  q  and  the  degree  of  decoherence  were  studied  in  microwave  billiards  where  a  non-zero  imaginary  part  of 
q  was  declared.  Moreover,  it  was  suggested  that  the  trajectory  of  q  in  its  complex  plane  could  be  used  to 
probe  whether  decoherence  is  caused  by  dissipation  or  dephasing.  Complex  q  parameter  was  used  in  other 
contexts  as  well,  such  as  the  second-order  effects  in  helium  auto-ionization  excited  by  electron  impact  and 
ultrashort  laser  excitation  in  the  bismuth  single  crystal. 

During  the  project  period,  we  established  the  universality  of  the  asymmetric  Fano  resonance  profile  and 
addressed  the  question  of  whether  the  fundamental  q  parameter  should  be  complex  or  real.  To  be  concrete, 
we  focused  on  low-dimensional  quantum-transport  systems  exemplified  by  quantum  dots  or  quantum  point 
contacts.  Despite  the  relatively  long  history  of  research  on  Fano  resonance,  for  quantum  transport  systems 
the  asymmetric  profile  had  been  sfudied  much  later.  Our  approach  to  probing  into  the  nature  of  the  q  parame¬ 
ter  consisted  of  two  steps.  First,  by  using  the  non-equilibrium  Green’s  function  to  calculate,  for  all  scattering 
channels,  the  transmission  as  a  function  of  the  Fermi  energy,  we  derived  a  formula  for  the  resonance  profile 
and  verified  if  numerically  using  graphene  quantum  dot  systems.  The  key  approximation  employed  is  that 
the  self-energy  terms  and  the  coupling  functions  are  slow  variables,  so  for  energy  near  an  isolated  reso¬ 
nance,  the  Green’s  function  can  be  decomposed  into  a  fast  and  a  slowly  varying  components.  The  formula 
also  leads  to  an  expression  of  the  width  of  the  resonance.  Second,  we  showed  that  our  formula  bridges 
naturally  the  conventional  Fano  formulas  with  complex  and  real  q  parameters.  In  particular,  we  found  a 
simple  mathematical  transform  that  can  convert  the  Fano  formula  with  complex  q  parameter  to  our  formula, 
and  another  transform  that  turns  our  formula  into  the  Fano  formula  with  real  q  parameter.  The  implication 
is  that,  for  any  experimental  situation  the  conventional  Fano  formulas  with  complex  or  real  q  parameter  are 
equally  applicable,  given  that  an  offset  can  be  added  to  the  real  q  Fano  formula.  In  our  analysis,  no  detailed 
information  about  the  specific  system  is  required.  Thus,  our  formula  and  its  direct  consequence  hold  for  any 
coherent  transport  dynamics  for  both  bosons  and  fermions.  Especially,  for  bosons,  i.e.,  phonons,  the  con¬ 
trol  parameter  is  not  the  energy  but  the  frequency,  thus  requiring  only  a  straightforward  modification  in  our 
formula.  For  fermion  transport,  our  formula  is  valid  for  small-seale  electronic  devices  such  as  quantum-dot 
systems,  quantum  point  contacts,  nanoscale  heterostructures,  and  single-molecule  transport  devices. 

Details  of  this  work  can  be  found  in 

•  L.  Huang,  Y.-C.  Lai,  H.-G.  Luo,  and  C.  Grebogi,  “Universal  formalism  of  Lano  resonance,”  AIP 
Advances  5,  017137,  1-18  (2015). 


3.3  Optical  branched  waves,  optomechanical  systems,  and  fluid  turbulence 

3.3.1  Emergence  of  scaling  associated  with  complex  branched  wave  structures  in  optical  medium 

When  waves  propagate  through  random  media,  extreme  events  and  complex  structures  such  as  rogue  waves 
and  branched,  fractal -like  wave  patterns  can  form.  There  is  a  great  interest  in  complex  wave  phenomena  due 
to  their  occurrences  in  a  host  of  physical  systems.  Lor  example,  in  oceanography,  rogue  waves  are  an  issue 
of  great  concern.  In  the  past  fifteen  years  or  so  there  had  been  experimental  and  theoretical  studies  of  rogue 
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Figure  10:  Fractal-like  branched  wave 
patterns  in  optical  media.  For  a  rect¬ 
angular  optical  medium  of  size  35/jm  x 
lQfm\  with  N  =  300  spatially  localized, 
Gaussian-shaped  scatterers,  typical  spatial 
distributions  of  the  magnetic  field  strength 
|H|  of  the  scattered  waves.  In  each  case,  an 
external  polarized,  monochromatic,  Gaus¬ 
sian  wave  of  width  w  =  ljum  and  unit 
intensity  is  sent  from  the  top  of  the  re¬ 
gion.  The  spatial  distribution  function  of 
the  refractive  index  for  the  whole  system  is 
n(r)  =  nQ  +  Zf=i  exp[-  |r-  r, 12/(20^)], 
where  An,  is  the  magnitude  of  the  refrac¬ 
tive  index  of  the  ith  scatterer  relative  to 
that  of  the  medium,  and  r  =  (x,  j)  is  a  two- 
dimensional  vector.  Each  scatterer  is  char¬ 
acterized  by  a  Gaussian-shaped  refractive 
index  profile,  whose  effective  radius  is  c. 
For  panels  (a)  and  (b),  the  refractive-index 
variations  are  negative:  An,-  =  —0.5.  The 
variation  is  positive  for  panel  (c):  An,-  = 
0.5.  For  panel  (d),  An,  is  randomly  selected 
from  the  range  [—0.5, 0.5].  Other  param¬ 
eters  are  the  same  for  all  panels:  no  =  1, 
X  =  1/rm,  and  G  =  0.22. 


waves  arising  from  long-range  acoustic  wave  propagation  through  ocean’s  sound  channel,  as  well  as  large- 
scale  experiments  on  directional  ocean  waves  to  probe  the  physical  and  dynamical  origin  of  these  extreme 
waves.  Extreme  events  and  complex  wave  patterns  were  also  identified  in  many  other  physical  situations 
such  as  light  propagation  in  doped  fibers,  acousfic  turbulence  in  superfluid  helium,  resonances  in  nonlinear 
optical  cavities,  linear  light-wave  propagation  in  multi-mode  glass  fiber,  and  elecfronic  fransport  in  semi- 
conducfor  two-dimensional  electron  gas  (2DEG)  systems.  Despite  previous  efforts,  an  accepted,  relatively 
complete  understanding  of  complex  extreme  waves  at  the  level  of  fundamental  physics  was  lacking. 

To  illustrate  the  extent  to  which  complex  branched  wave  patterns  were  previously  understood,  we  choose 
electronic  transport  in  2DEG  systems  as  an  example.  It  was  found  by  a  Flarvard  group  in  2001  that  electron 
flows  from  a  quanfum  point  contact  exhibit  a  branched  or  fractal-like  behavior  with  highly  non-uniform 
amplitude  distribution  in  the  physical  space.  The  observed  separated,  narrow  strands  of  greatly  enhanced 
electron  wave  intensities  were  argued  to  be  caused  by  random  background  potentials  and  quantum  coherent 
phase  interference  among  the  electronic  wave  functions.  Subsequently  a  theory  was  proposed  to  predict 
the  statistical  distribution  of  the  intensities  of  branched  electron  flows  in  fhe  presence  of  weak,  correlated 
Gaussian  random  potentials. 

The  generic  origin  of  wave  branching  behavior  has  been  a  matter  of  active  debate.  A  tacit  assumption 
in  most  previous  investigations  was  nonlinearity  in  the  underlying  medium.  In  particular,  it  had  been  be¬ 
lieved  that  the  existence  of  many  uncorrelated,  spatially  randomly  distributed  wave  elements  is  key  to  the 
occurrence  of  these  exotic  wave  patterns.  These  elements  can  be,  for  example,  solitons  in  nonlinear  systems. 
However,  it  was  demonstrated  experimentally  in  a  microwave  system  and  in  a  multi-mode  optical  fiber  fhaf 
branched  wave  patterns  can  occur  even  in  fhe  absence  of  nonlinearity.  In  fact,  in  the  latter  case,  granularity 
of  light  speckles  at  the  fiber  exif  and  inhomogeneify  in  the  spatial  clustering  of  the  speckle  patterns  are 
speculated  to  be  the  two  ingredients  that  trigger  complex  wave  patterns.  These  works  thus  demonstrated 
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Figure  1 1 :  Intensity  scaling  of  fractal¬ 
like,  branched  wave  patterns  in  opti¬ 
cal  media,  (a)  Distribution  of  50  ran¬ 
dom  scatterers  in  the  35fjm  x  70/rm  uni¬ 
form  medium,  (b)  Branched  wave  mag¬ 
nitude  from  FDFD  simulation.  The  inci¬ 
dent  wave  is  uniform  with  width  X  and  sent 
in  the  -by  direction,  where  X  =  l/7m,no  = 
l,An  =  —0.5, a  =  0.23/im.  (c)  Numeri¬ 
cally  obtained  scaling  law  of  the  intensity 
distribution  (blue  circles)  on  a  logarithmic 
scale  and  the  theoretical  prediction  (solid 
line  ^  /^^).  There  is  an  exponential  cut¬ 
off  in  the  high-intensity  regime,  as  shown 
in  the  inset. 


that  nonlinearity  is  not  absolutely  essential  for  the  emergence  of  these  extreme  waves.  An  issue  of  signifi¬ 
cant  theoretical  and  experimental  interest  concerns  thus  a  minimal,  physical  model  that  can  generate  robust 
branched  wave  patterns,  so  that  their  generic  and  physical  origin  may  be  elucidated.  A  related  issue  is  the 
statistical  properties  of  these  waves.  In  this  regard,  a  general  observation  in  all  contexts  where  branched 
wave  structures  arise  is  the  non-Gaussian  statistics  of  the  wave  amplitude.  Typically  there  is  a  long  tail 
in  the  probability  density  function,  which  characterizes  the  extreme  intensity  of  the  waves.  An  essential 
requirement  for  a  valid  minimal  model  of  branched  wave  patterns  is  thus  that  it  generate  the  universally 
observed  long-tail  distribution  in  the  wave  intensity. 

During  the  project  period,  we  proposed  a  class  of  minimal  models  for  branched  wave  patterns  in  the 
context  of  wave  propagation  in  two-dimensional  optical  medium.  The  model  contains  two  basic  physical 
elements:  (1)  a  uniform  medium  of  finite  size  and  (2)  spatially  localized  scatterers  randomly  distributed  in 
the  medium,  the  refractive  indices  of  which  deviate  from  that  of  the  background  medium.  The  deviations  can 
occur  in  both  ways  which,  in  the  case  of  negative  deviation,  may  correspond  to  scatterers  that  are  effectively 
negative-indexed,  or  metamaterials.  The  second  element  is  required  for  generating  dynamics  beyond  simple 
linear  wave  propagation.  We  demonstrated  that  such  a  minimal  model  can  generate  robust  branched  wave 
patterns,  regardless  of  the  detailed  distribution  of  the  refractive-index  deviations  associated  with  the  random 
scatterers.  We  then  developed  a  comprehensive  analytic  treatment  of  the  minimal  model  by  focusing  on  the 
theoretical  derivation  of  the  power-law  type  of  long-tail  type  of  distribution  in  the  wave  intensities. 

More  specifically,  we  considered  the  setting  where  a  polarized  monochromatic  light  propagates  in  a 
dielectric  optical  medium  with  random  structural  imperfections  characterized  by  random  refractive-index 
disorders  (scatterers)  of  size  comparable  to  the  wavelength.  We  then  employed  the  standard  finite-difference 
frequency-domain  (FDFD)  method  to  calculate  the  intensity  of  the  scattered  field  through  multiple  scatterers. 


18 


In  the  weak  scattering  limit,  i.e.,  when  the  wavelength  X  is  much  smaller  than  the  mean  free  path  /,  we 
obtained  striking  branching  flow  structures  of  propagating  light,  similar  to  those  observed  in  the  2DEG  and 
microwave  transport  experiments.  Examples  of  branched  wave  patterns  are  shown  in  Eig.  10.  As  the  spatial 
density  of  the  scatterers  is  increased,  the  intensity  patterns  exhibit  more  pronounced  fractal-like  behavior, 
where  branches  of  extraordinarily  high  intensities  tend  to  enhance  themselves  when  forking  into  narrower 
and  even  smaller  paths.  Anderson  localization  of  light  is  also  observed  as  the  mean  free  path  approaches  the 
strong  scattering  limit  (/  ~  X).  Our  extensive  numerical  computation  confirmed  that  the  branched  structure 
can  result  from  the  caustics  of  the  flow  rather  than  the  valleys  of  the  random  scatterers.  We  found  that 
branched  waves  generically  arise  in  the  regime  between  weak  scattering  and  strong  localization  of  light 
waves. 

In  order  to  obtain  a  comprehensive  understanding  of  the  occurrence  of  branched  waves  in  optical  media 
and  also  to  uncover  their  statistics,  we  developed  a  detailed  analytic  theory.  Utilizing  the  Green’s  function 
method,  we  treated  the  scattering  of  two-dimensional  polarized  light  wave  off  a  single  scatterer  and  obtained 
a  theoretical  explanation  for  the  reason  why  the  wavelength  needs  to  be  comparable  to  the  size  of  the  scatterer 
in  order  for  noticeable  large  fluctuations  of  branching  strands  to  be  observed.  In  contrast  to  the  existing 
theory  which  deals  with  stretches  and  folds  in  classical  ray  dynamics  in  a  concrete  and  somewhat  abstract 
manner,  our  theory  enables  us  to  visualize  the  branching  flow  structures  of  the  scattered  light  intensities  in 
different  angular  directions.  Based  on  the  results  from  scattering  off  a  single  scatterer,  we  next  extended 
our  treatment  to  multiple  scatterers.  In  this  case,  coherent  backscattering  and  recurrent  multiple  scattering 
become  important,  and  they  together  contribute  dominantly  to  the  formation  of  extremely  large  amplitude 
events.  This  unstable  branch  stretching  and  accumulation  process  is  highly  sensitive  to  the  scatterers’  spatial 
distribution  and  thus  is  critical  to  the  formation  of  fractal-like  wave  patterns.  Because  of  the  large  intensity 
fluctuations  caused  by  wave  interference  and  the  complexity  of  random- scatterer  configuration,  it  is  essential 
to  focus  on  the  statistical  distribution  function  of  light  intensities.  We  succeeded  in  deriving  a  formula  for 
the  distribution  function  of  high  intensities,  which  follows  an  algebraic  scaling  law  in  the  weak-scattering 
limit.  An  example  of  the  algebraic  intensity  distribution  is  shown  in  Eig.  11.  The  physical  significance  of 
the  algebraic  scaling  means  that,  associated  with  the  branched  waves,  there  are  points  in  the  space  at  which 
exceedingly  large  intensities  can  arise,  the  “hot”  spots,  in  contrast  to  situations  governed  by  Gaussian  type 
of  intensity  distributions. 

Details  of  this  work  can  be  found  in 

•  X.  Ni,  Y.-C.  Lai,  and  W.-X.  Wang,  “Emergence  of  scaling  associated  with  complex  branched  wave 
structures  in  optical  medium,”  Chaos  22,  043116,  1-11  (2012). 


3.3.2  Enhancement  of  quantum  entanglement  by  nonlinear  dynamics  in  optomechanical  systems 

When  a  quantum  system  contains  several  degrees  of  freedom,  entanglement  among  them  can  arise.  Quantum 
entanglement  is  a  form  of  quantum  superposition,  which  can  be  best  understood  by  considering  a  pair  of 
particles  such  as  electrons.  When  two  such  particles  are  entangled,  each  particle  cannot  be  fully  described 
without  considering  the  other  because  they  both  belong  to  a  single  quantum  superposition  state.  Without 
external  disturbance,  e.g.,  a  measurement,  the  particles  remain  in  the  entangled  state.  In  order  to  determine 
the  value  of  some  physical  quantity  associated  with  one  particle  in  the  entangled  state,  it  is  necessary  to 
apply  perturbation  to  the  system,  after  which  the  state  of  the  particle  takes  on  a  definite  value.  However,  the 
other  particle  in  the  originally  entangled  pair  will  have  the  corresponding  correlated  value  for  any  subsequent 
time,  regardless  of  the  physical  distance  between  the  two  particles.  The  remarkable  phenomenon  of  quantum 
entanglement  has  been  observed  in  experiments  of  large  molecules  and  even  small  diamonds,  and  it  is  the 
foundation  of  modern  quantum  computing  and  information  science. 
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Figure  12:  Nonlinear  dynamics  in  op¬ 
tomechanical  systems,  (a)  Schematic  il¬ 
lustration  of  a  typical  optomechanical  sys¬ 
tem  driven  by  a  single-mode  laser.  The 
laser  beam  enters  the  optical  cavity  through 
a  fixed  mirror  but  photons  inside  the  cav¬ 
ity  can  also  decay  through  the  mirror.  The 
movable  mirror  is  totally  reflecting  and 
linked  to  a  mechanical  oscillator.  (b,c) 
Time  series  of  the  dimensionless  canonical 
coordinate  q  for  situations  where  the  clas¬ 
sical  dynamics  are  periodic  and  quasiperi- 
odic,  respectively. 


While  the  principle  of  superposition  is  defined  with  respect  to  linear  systems,  our  physical  world,  when 
viewed  classically,  is  nonlinear.  A  nonlinear  dynamical  system  can  exhibit  all  kinds  of  interesting  phe¬ 
nomena  such  as  periodic  oscillations,  quasiperiodicity  and  chaos,  and  they  can  have  distinct  fingerprints 
or  manifestations  when  the  system  is  treated  quantum-mechanically.  Especially,  the  studies  of  quantum 
manifestations  of  classical  chaos  constitute  the  field  of  quantum  chaos,  and  there  have  been  tremendous 
efforts  in  the  past  four  decades  in  this  field.  However,  there  is  no  classical  correspondence  of  quantum  en¬ 
tanglement  in  the  strict  sense,  as  the  measurement  of  a  physical  state  at  one  place  immediately  influences 
the  measurement  at  the  other.  In  view  of  the  ubiquity  of  nonlinear  dynamics  in  physical  systems  and  the 
fundamental  importance  of  quantum  entanglement,  curiosity  demands  that  we  ask  the  following  question: 
What  is  the  interplay  between  nonlinear  dynamics  and  quantum  entanglement?  Or  more  explicitly,  when  the 
corresponding  classical  system  is  nonlinear,  what  are  the  signatures  on  quantum  entanglement?  This  issue 
may  have  significant  practical  implications  in  developing  devices  and  systems  for  quantum  computing  and 
quantum  information  processing. 

During  the  project  period,  we  addressed  the  nonlinear-dynamics/quantum-entanglement  issue  by  using 
optomechanical  systems,  a  field  of  intense  recent  investigation.  As  shown  in  Fig.  12,  an  optomechanical 
system  consists  of  an  optical  cavity  and  a  nanoscale  mechanical  oscillator,  such  as  a  cantilever.  When 
a  laser  beam  is  introduced  into  the  cavity,  a  resonant  optical  field  is  established  that  exerts  a  radiation 
force  on  the  mechanical  cantilever,  causing  it  to  oscillate.  The  mechanical  oscillations  in  turn  modulate  the 
length  of  the  optical  cavity,  hence  its  resonant  frequency.  There  is  thus  coupling  between  the  optical  and 
the  mechanical  degrees  of  freedom.  This  coupling,  or  interaction,  can  lead  to  cooling  of  the  mechanical 
oscillator  toward  the  quantum  ground  state,  a  topic  of  great  recent  interest.  The  optomechanical  coupling 
thus  provides  a  straightforward  way  to  entangle  the  optical  with  the  mechanical  modes.  Previous  works  on 
quantum  entanglement  in  optomechanical  systems  focused  on  the  situations  where  the  classical  dynamics 
are  either  steady-state  or  periodic  oscillations.  Our  goal  was  to  search  for  nonlinear  dynamical  behaviors 
beyond  steady-state  solutions  and  periodic  oscillations,  and  to  investigate  the  effects  of  such  behaviors  on 
quantum  entanglement. 

Our  main  findings  were  the  following.  In  an  experimentally  realizable  parameter  regime,  as  the  power 
of  the  driving  laser  is  increased,  there  is  a  transition  from  periodic  to  quasiperiodic  motions,  where  in  the 
latter,  the  system  is  strongly  nonlinear  with  two  incommensurate  frequencies.  Entanglement  is  enhanced 
towards  the  transition,  vanishes  as  the  transition  point  is  reached,  but  is  restored  abruptly  after  the  transition 
and  continues  to  be  enhanced  as  the  system  evolves  deeply  into  the  quasiperiodic  regime.  Representative 
results  are  shown  in  Fig.  13.  A  surprising  result  is  that,  with  respect  to  time  evolution,  there  is  a  direct  corre¬ 
spondence  of  quantum  entanglement  to  classical  nonlinear  dynamics.  In  particular,  for  classically  periodic 
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Figure  13:  Enhancement  of  qnantnm 
entanglement  by  nonlinear  dynamics  in 
optomechanical  systems,  (a)  Typical  bi¬ 
furcation  diagram  of  the  generic  optome¬ 
chanical  system  in  Fig.  12(a),  where  the 
asymptotic  values  of  the  dimensionless 
canonical  coordinate  q  are  plotted  versus 
the  normalized  external  laser  power  P.  The 
purple  dashed  line  indicates  the  threshold 
driving  power,  at  which  a  transition  from 
periodic  to  quasiperiodic  motions  occur. 
The  cavity  length  is  fo  =  25mm  with  finesse 
F  =  6  X  10^.  (b)  The  entire  spectrum  of 
the  Lyapunov  exponents  versus  P,  and  (c) 
the  maximum  of  the  pseudo-entanglement 
measure  Ep  versus  P. 


dynamics,  the  time  evolution  of  the  entanglement  measure  is  also  periodic,  but  when  the  classical  system 
enters  into  the  quasiperiodic  regime,  the  quantum  entanglement  measure  exhibits  a  beats-like  behavior  with 
two  distinct  frequencies.  The  entanglement,  especially  when  the  classical  system  is  quasiperiodic,  is  robust 
with  respect  to  temperature  variations.  Thus,  from  the  perspective  of  entanglement,  there  is  a  clear  quantum 
manifestation  and  signature  of  nonlinear  dynamics  in  optomechanical  systems. 

Details  of  this  work  can  be  found  in 

•  G.-L.  Wang,  L.  Fluang,  Y.-C.  Lai,  and  C.  Grebogi,  “Nonlinear  dynamics  and  quantum  entanglement 
in  optomechanical  systems,”  Physical  Review  Letters  112,  110406,  1-6  (2014). 


3.3.3  Low  Reynolds-number  turbulence  in  ferrofluids 

A  basic  understanding  of  fluid  turbulence  is  that  it  typically  occurs  at  high  values  of  the  Reynolds  number,  a 
dimensionless  quantity  defined  as  Re  =  lu/v,  where  /  and  u  are  the  characteristic  length  scale  and  a  typical 
velocity  of  the  flow,  respectively,  and  v  is  the  kinematic  viscosity.  Take,  for  example,  the  paradigmatic 
setting  of  a  uniform  flow  of  velocity  u  flowing  past  a  cylinder  of  diameter  1.  When  Re  is  of  the  order 
of  tens,  the  flow  is  regular.  For  Re  in  the  hundreds,  von  Karman  vortex  street  forms  behind  the  cylinder, 
breaking  certain  symmetries  of  the  system.  Fully  developed  turbulence,  in  which  the  broken  symmetries 
are  restored,  occurs  at  very  high  Reynolds  number,  typically  in  the  thousands.  Because  of  the  requirement 
of  high  Reynolds  numbers,  study  of  turbulence  has  historically  been  known  as  an  extremely  challenging 
problem,  both  experimentally  where  flow  systems  of  enormous  size  and/or  high  velocity  are  required  and 
computationally  where  unconventionally  high  resolution  in  the  numerical  integration  of  the  Navier-Stokes 
equation  is  needed.  It  is  thus  desirable  that  fluid  turbulence  can  emerge  in  physical  flows  of  relatively  low 
Reynolds  numbers.  At  a  fundamental  level,  such  a  flow  system  may  present  us  with  a  new  paradigm  that  can 
potentially  yield  significantly  new  insights  into  the  emergence,  characterization,  mechanism,  and  control  of 
fluid  turbulence. 

During  the  project  period,  we  investigated  the  Taylor-Couette  flow  in  finite  systems  (e.g.,  aspect  ratio 
F  =  20)  where  a  rotating  ferrofluid  is  confined  by  axial  end  walls,  i.e.,  non-rotating  lids,  in  the  presence  of 
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(d)  turbulence 
s  =  0.9 


Figure  14:  Transition  to  turbulence  in 
a  rotating  ferrofluid.  (a-d)  For  four  val¬ 
ues  of  the  magnetic  parameter  (Niklas  pa¬ 
rameter)  Sx  corresponding  to  WVF2  (wavy 
vortices  with  a  two-fold  symmetry),  WVF; 
(wavy  vortices  with  a  tilting  pattern), 
SWOp  (standing  waves  with  axial  oscil¬ 
lations),  and  turbulence  regimes,  respec¬ 
tively,  isosurfaces  of  azimuthal  vorticity 
and  contours  of  the  radial  velocity  on  an 
unrolled  cylindrical  surface  in  the  annulus 
at  midgap.  Red  (dark  gray)  and  yellow 
(light  gray)  colors  denote  isosurfaces,  and 
inflow  and  outflow  for  contour  plots,  re¬ 
spectively.  The  patterns  in  (c,d)  are  snap¬ 
shots  due  to  time  dependence  of  the  corre¬ 
sponding  flow. 


an  external  magnetic  field.  The  classic  Taylor-Couette  flow  of  non-ferrofluid  has  been  a  computational  and 
experimental  paradigm  to  investigate  a  variety  of  nonlinear  and  complex  dynamical  phenomena,  including 
transition  to  turbulence  at  high  Reynolds  numbers.  The  corresponding  ferrofluid  system  we  sfudied  con¬ 
sists  of  two  independently  rotating,  concentric  cylinders  with  viscous  ferrofluid  filled  in  befween,  which  has 
embedded  within  itself  artificially  dispersed  magnetized  nanoparticles.  In  absence  of  the  external  magnetic 
field,  the  magnetic  moments  of  the  nanoparticles  are  randomly  oriented,  leading  to  zero  net  magnetization 
for  the  entire  fluid.  In  this  case,  the  magnetized  nanoparticles  have  little  effect  on  the  physical  proper¬ 
ties  of  the  fluid  such  as  density  and  viscosity.  However,  when  a  transverse  magnetic  field  is  applied,  the 
physical  properties  of  the  fluid  are  significanfly  modified,  leading  to  drastic  changes  in  the  underlying  hy¬ 
drodynamics.  For  example,  for  systems  of  rotating  ferrofluid,  an  external  magnetic  field  can  stabilize  regular 
dynamical  states  and  induce  dramatic  changes  in  the  flow  topology.  In  general,  the  magnetic  field  can  be 
used  effeclively  as  a  control  or  bifurcation  parameter  of  the  system,  whose  change  can  lead  to  character¬ 
istically  distinct  types  of  hydrodynamical  behaviors.  Existing  works  on  rotating  ferrofluids  however,  were 
mostly  concerned  with  steady  time -independent  flows.  Time-dependent  ferrofluid  flows  were  investigated 
only  recently  but  in  the  non-turbulent  regime. 

The  interaction  between  ferrofluid  and  magnetic  field  leads  to  additional  terms  in  the  Navier-Stokes 
equation.  Our  extensive  computations  revealed  a  sequence  of  bifurcations  leading  to  time-dependent  flow 
solutions  such  as  standing  waves  with  periodic  or  quasiperiodic  oscillations,  and  turbulence.  Surprisingly, 
we  found  that  turbulence  can  occur  for  Reynolds  numbers  at  least  one  order  of  magnitude  smaller  than  those 
required  for  turbulence  to  arise  in  conventional  fluids.  The  occurrence  of  turbulence  was  ascertained  by  a 
bifurcation  analysis  and  by  examining  the  characteristics  of  the  physical  quantities  such  as  the  energy,  the 
wave  number,  and  the  angular  momentum.  We  also  found  that  the  onset  of  turbulence  can  be  determined 
accurately,  in  contrast  to  classical  fluid  turbulence  where  such  a  determination  is  typically  qualitative  and 
involves  a  high  degree  of  uncertainty.  Our  findings  have  fhe  following  implicafions:  (1)  ferrofluids  under 
magnetic  field  is  a  new  paradigm  for  investigating  turbulence,  especially  experimentally  where  the  study 
can  be  greatly  facilitated  due  to  the  dramatic  relaxation  in  the  Reynolds-number  requirement,  (2)  the  critical 
magnetic-field  strength  for  the  onset  of  turbulence  can  be  pinned  down  precisely,  possibly  leading  to  deeper 
insights  into  the  physical  and  dynamical  origins  of  the  transition,  and  (3)  turbulence  can  be  controlled 
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externally,  e.g.,  by  an  external  magnetie  field. 

Details  of  this  work  can  be  found  in 

•  S.  Altmeyer,  Y.-H.  Do,  and  Y.-C.  Lai,  “Transition  to  turbulence  in  Taylor-Couette  ferrofluidic  flow,” 
Nature  Scientific  Reports  5,  10781,  1-12  (2015). 

3.4  Superpersistent  currents  and  whispering  gallery  modes  in  relativistic  quantum  chaotic 
systems 

A  remarkable  phenomenon  in  the  quantum  world  is  persistent  currents  (PCs),  permanent  currents  without 
any  external  source,  which  are  generated  by  the  Aharonov-Bohm  (AB)  effect  in  non- superconducting  sys¬ 
tems.  PCs  have  been  observed  experimentally  in  metallic  and  semiconductor  rings  in  the  mesoscopic  regime. 
Theoretical  efforts  have  been  focused  on  the  effects  of  bulk  disorders,  electron-electron  interactions,  spin- 
orbital  interactions,  and  electromagnetic  radiation  on  PCs,  typically  studied  in  the  diffusive  regime  using 
idealized  circular-symmetric  rings  and  cylinders.  Rapid  advances  in  nanotechnology  have  made  it  feasible 
to  fabricate  mesoscopic  devices  with  mean  free  path  larger  than  their  sizes  at  sufficiently  low  temperatures 
(the  ballistic  transport  regime).  The  AB  system  can  thus  be  modeled  as  a  quantum  ballistic  billiard  in  which 
the  particles  are  scattered  at  the  boundaries  of  the  domain.  As  a  result,  the  boundary  shape  becomes  highly 
relevant.  In  experiments,  uncontrollable  boundary  imperfections  are  inevitable  even  when  there  are  no  bulk 
disorders.  It  is  thus  of  interest  to  study  the  effects  of  boundaries,  e.g.,  those  that  generate  chaos  in  the 
classical  limit,  on  PCs.  In  general,  an  asymmetric  boundary  destroys  angular  momentum  conservation  and 
introduces  irregular  scattering.  Theoretical  and  experimental  studies  have  shown  that,  similar  to  the  effects 
of  bulk  disorder,  symmetry  breaking  of  the  boundary  can  result  in  opening  of  gaps  at  the  degeneracy  points 
of  the  energy  levels,  leading  to  level  repulsion,  a  typical  manifestation  of  classical  chaos.  Energy  gap  open¬ 
ing  can  diminish  AB  oscillations  through  pinning  of  the  corresponding  states,  leading  to  vanishing  PCs. 
Since  fully  chaotic  domains  are  associated  with  a  strong  degree  of  symmetry  breaking,  PCs  are  not  expected 
to  arise.  In  nonrelativistic  quantum  systems  governed  by  the  Schrddinger  equation,  PCs  are  thm  fragile. 

We  addressed  the  question  of  whether,  in  relativistic  quantum  systems,  PCs  can  arise  and  sustain  in 
the  presence  of  symmetry-breaking  perturbations.  Besides  the  importance  of  this  question  to  fundamental 
physics,  our  work  was  motivated  by  the  tremendous  recent  research  on  two-dimensional  Dirac  materials  such 
as  graphene,  topological  insulators,  molybdenum  disulfide  (M0S2),  HITP  [Ni3(HITP)2],  and  topological 
Dirac  semimetals.  The  physics  of  these  materials  is  governed  by  the  Dirac  equation.  This  is  thus  interest  in 
investigating  relativistic  PCs  in  Dirac  fermion  systems. 

Existing  theoretical  works  on  relativistic  PCs,  however,  assumed  idealized  circular-symmetric  rings  in 
the  ballistic  limit.  Whether  AB  oscillations  and  consequently  relativistic  PCs  can  exist  in  asymmetric  rings 
that  exhibit  chaos  in  the  classical  limit  is  a  fundamental  question.  Our  finding  was  that,  even  in  the  presence 
of  significant  boundary  deformations  so  that  the  classical  dynamics  becomes  fully  chaotic,  robust  PCs  can 
occur  in  relativistic  quantum,  Dirac  fermion  systems,  henceforth  the  term  superpersistent  currents  (SPCs). 
A  more  striking  finding  is  that  SPCs  are  generated  by  localized  states  at  the  domain  boundaries,  which 
are  effectively  chaotic  Dirac  whispering  gallery  modes  (WGMs)  that  carry  larger  angular  momenta.  While 
WGMs  are  common  in  photonic  systems,  its  emergence  in  electronic  systems,  especially  in  relativistic 
quantum  systems,  is  rare  and  surprising.  We  developed  a  physical  understanding  of  the  counterintuitive 
phenomenon  of  SPCs  by  analytically  exploiting  the  properties  of  the  spinor  wavefunctions  in  an  idealized 
relativistic  quantum  system. 

The  significance  of  our  results  lies  in  the  perspective  of  observing  SPCs  in  the  presence  of  strong  ran¬ 
dom  scattering,  implying  that  they  may  occur  in  systems  of  size  beyond  the  mesoscopic  limit.  This  can 
potentially  be  a  relativistic  quantum  phenomenon  occurring  at  relatively  large  scales.  There  can  be  signif¬ 
icant  applications  of  SPCs  in  quantum  information  processing,  for  which  we  presented  a  scheme  of  Dirac 
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WGM-based  qubit.  We  note  that,  previous  experimental  studies  of  AB  oscillations  in  graphene  and  topo¬ 
logical  insulators  made  it  possible  to  experimentally  test  the  phenomena  of  Dirac  WGMs  and  SPCs  that  we 
predicted.  To  motivate  experimental  verihcation,  we  proposed  a  feasible  scheme  using  three-dimensional 
topological  insulators. 

Contrast  of  eigen-wavefunctions  for  the  Schrodinger  particle  and  Dirac  fermion  is  demonstrated  in 
Fig.  15.  For  low  energy  levels,  the  Schrodinger  particle  is  strongly  localized  throughout  the  domain,  as 
shown  in  Figs.  15(a-c),  leading  to  a  flat  energy-flux  dispersion.  However,  the  Dirac  fermion  typically  travels 
around  the  ring’s  boundaries,  forming  relativistic  WGMs  that  persist  under  irregular  boundary  scattering  due 
to  chaos  and  are  magnetic-flux  dependent,  as  shown  in  Fig.  15(d-f).  Conventional  wisdom  for  Schrodinger 
particle  stipulates  that  asymmetry  in  the  domain  geometry  can  mix/couple  well-defined  angular  momentum 
states,  opening  energy  gaps  and  leading  to  localization  of  lower  states  in  the  entire  domain  region,  so  AB  os¬ 
cillations  would  vanish,  as  demonstrated  both  theoretically  and  experimentally.  However,  for  Dirac  fermion, 
this  picture  breaks  down  -  there  are  robust  AB  oscillations  even  in  the  fully  chaotic  domain  and  the  particle 
tends  to  execute  motions  corresponding  to  WGMs,  leading  to  SPCs. 

Figure  15:  Superpersistent  current  in 
two-dimensional  Dirac  fermion  systems. 

Probability  distribution  of  the  10th  eigen¬ 
state  for  (a-c)  nonrelativistic  and  (d-f)  rel¬ 
ativistic  AB  chaotic  billiard  for  magnetic- 
flux  parameter  a  =  —1/4,0, 1/4,  respec¬ 
tively. 


•  H.-Y.  Xu,  L.  Huang,  Y.-C.  Lai,  and  C.  Grebogi,  “Superpersistent  currents  and  whispering  gallery 
modes  in  relativistic  quantum  chaotic  systems,”  Nature  Scientific  Reports  5,  8963,  1-8  (2015). 
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